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INTERACTION BETWEEN A SUPERSONIC STREAM AND A PARALLEL 
SUBSONIC STREAM BOUNDED BY FLUID AT REST 
By Herbert S. Ribner and E. Leonard Amoff 


SUMMARY 

In a simplified inviscid model of shock-wave boundary-layer inter- 
action, Tsien and Finston have replaced the boundary layer by a unif orm 
subsonic stream, bounded on one side by a solid wall nnr) on the other 
side by the interface with a uniform supersonic stream of s emi - infinit e 
extent. Among other things, this model fails to simulate the separated 
region or "dead -air" bubble that generally appears in a laminar boundary 
layer subjected to an oblique incident shock wave of moderate strength. 

In order to introduce a main feature of such a dead-air region, the model 
has been modified herein by replacing the solid wall by an interface with 
fluid at rest. 

The presence of the boundary layer sandwiched between the outer 
supersonic flow' and the dead-air region Is found scarcely to modify the 
shape, in the vicinity of the shock, of the expansive "corner" turn that 
would exist if" the shock were incident directly on the dead-air region 
without the intermediary of the boundary layer; there are local distor- 
tions top and bottom, but these are reduced to negligible amounts several 
boundary-layer thicknesses to the left or right of the effective corner. 

In support of a phase of, the work of Lester Lees, it is concluded 
that in a more accurate treatment of the complete region of shock 
boundary-layer interaction the Prandtl boundary-layer equations nay be 
applied to the entire extent of the disturbed boundary layer, applying 
as a boundary condition a sudden turn of the displacement surface through 
twice the shock deflection angle at the point of shock incidence. There- 
by the flow details in the immediate vicinity of the shock will be some- 
what in error, but the over-all features of the interaction are capable 
of being given correctly. Present unknown elements in such an application 
appear to be the point of transition from lamimr to turbulent flow anfl 
the form of the turbulent equations where a separated bubble exists. 
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INTRODUCTION 

The complex problem of the interaction between shock waves and 
boundary layers in supersonic flow is not yet clearly understood, 
despite a number of theoretical studies (references 1 to 7). Among 
these, the analysis of Tsien and Finston (reference 2) is singled 
out here for further examination. In the Tsien-Finston model the 
boundary layer is simulated by a inaiform subsonic stream of finite 
width bounded on one side by a solid wall and on the other side by 
the interface with a uniform supersonic stream of semi-infinite 
extent. The fluid is assumed to be nonviscous and nonheat -conducting 
and the disturbances are assumed small. 

The boundary-layer perturbation calculated in reference 2 from this 
model has the general character shown in sketch 1. The experimental 
perturbation for a laminar boundary layer (e.g., references 8 and 9), 
on the other hand, has the general character shown in sketch 2. The 
resemblance is not close. A striking experimental feature is the pres- 
ence of a separated region, which appears for all but the weakest inci- 
dent shock strengths (references 3 and 9). -The separated region is not 
simulated in the Tsien-Finston model of shock boundary-layer interaction, 
and this constitutes an important shortcom in g of the model for laminar 
boundary layers. For turbulent boundary layers subjected to moderate 
shock strengths, on the other hand, there will be no separation, and 
the model may be quite good. 

The air in a separated bubble is substantially at res^ . An oblique 
shock incident directly on a dead -air region, without the intermediary 
of a boundary layer, is known to reflect as an expansion wave (see 
sketch 3) . Note that the interface in sketch 3 is deflected much like 
the potential flow just above the actual boundary layer (sketch 2), a 
feature completely unduplicated in the Tsien-Finston model (sketch l) . 

A natural inference is that the separated region dominates the behavior 
in the immediate vicinity of the shock, and the upper surface behaves, 
in that vicinity, substantially as if the entire boundary layer were 
replaced by a dead-air region. 

This tempting inference has been made to provide a boundary condi- 
tion in the work of Lees (reference 6) and of Ritter (reference 7), along 
with other simplifying assumptions; it is the purpose of the present 
paper to test its validity on a simplified model applicable only to the 
Immediate vicinity of the shock. For this purpose a modification of 
the Tsien-Finston model will serve: the solid lower wall of sketch 1 is 

replaced by a semi-infinite region of air at rest, or "dead air". Thus 
the boundary condition of zero vertical perturbation velocity at the 
wall is replaced by a condition of constant pressure along the boundary 
of the dead-air region; in all other respects there is no change from 
the original model. Experimentally (reference 9), an approx imat ion to 
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the condition of constant pressure is found over a substantial portion 
of the separated air bubble including the part directly under the inci- 
dent shockj the ultimate sharp pressure rise generally occurs somewhat 
downstream of the shock. 

It has been implicit in the discussion that attention was limited 
to the boundary layer adjacent to a solid body. Also of interest;, how- 
ever, is the boundary layer between a two-dimensional supersonic jet and 
the stationary fluid bounding it. The present modified flow model will 
describe, in idealized fashion, the interaction of an oblique shock wave 
(such as may emanate from a nozzle lip) with such a jet boundary layer. 

The model is, in fact, better adapted to such a fluid-to -fluid boundary 
layer than it ig to the fluid-to -wall boundary layer. 

The mathematical analysis herein of the modified model is necessarily 
parallel in a number of respects to the original work of Tsien and Finston. 
Nevertheless, it has been thought worthwhile to provide a complete account 
In an appendix for purposes of unity and clarity. Hie symbols used in this 
analysis are presented in appendix Aj the analysis, in appendix B. An 
abbreviated version of the analysis is given in the main text. The inves- 
tigation was conducted at the NACA Lewis laboratory. 


SYNOPSIS OF ANALYSIS 


Flow regions and governing equations . - In figure 1 are shown the 
supersonic region (region l] at the top, the- subsonic region (region 2) 
in the middle, and the dead-air region (region 3) at the bottom. A weak 
shock wave, or Mach wave, is incident on the interface between, regions 1 
and 2. Hie governing equations for the perturbation potential are 


Region 1: 



d 2 <fh d 2 ^ 

+ 1 ? = °' 




/ 2\ d 2 <Po d 2 ?2 

Region 2: ll - Mg ] g + g = °> Mg < 1 

\ ) Sx <5y 

In region 1 the perturbation potential includes a part (incoming wave) 
which represents directly the incident Mach wave and a part (outgoing 
waves) which arises from the consequent warping of the interface between 
regions 1 and 2. 


Boundary conditions at each interface . - The boundary conditions are: 
(a) The perturbation pressure must be the same on both sides of an inter- 
face (taken to be zero at the lower interface), and (b) the flow inclina- 
tion must be the same on both sides. In linearized theory, boundary con- 
dition (a) is essentially a condition on the axial perturbation velocity 
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u, and “bo undar y condition (b) is essentially a condition on the trans- 
verse perturbation velocity v. At the lower interface, condition (h) 
yields no information, since no direction can he assigned to the zero 
velocity in the bottom region; thus (a) and (b), as applied to both 
interfaces, yield three boundary conditions. 

Solution in general form. - The general solution for each region is 
expressed in the form of a Fourier integral; such a form is capable of 
satisfying arbitrary boundary conditions . The variable of integration 
is X, and terms like A-^X) cos Xx, Bp(X) sin Xx, and so forth, appear; 

there are six Fourier coefficients Ap(X), Bj_(X), EgCX), G^CX), 

and HgCX) to be determined . 

Application of boundary conditions . - The velocity components in 
Fourier integral form are now substituted into the three boundary con- 
ditions; equating correspond in g coefficients of sine terms end cosine 
terms results, in effect, in six equations. Two of the equations yield 
E 2 (X) = F 2 (X) = 0; thus there remain four equations which are solved 
simultaneously for the four remaining Fourier coefficients A^(X) , B-^(X), 
G 2 (X), and H 2 (X). 

Solution for pressure along upper interface . - The perturbation 
pressure along the upper interface is proportional to u 2 there; the 
expression is of the form of a Fourier integral involving G 2 (X) and 
H 2 (X) , which are now known. Evaluation and simplification of the 
integral is complicated and tedious; the final result is expressed in 
series form and evaluated numerically. In this evaluation, functions 
tabulated by Tsien and Finston are employed. 

Solution for displacement of each interface . - The local displace- 
ment t)(x) of an interface from its undisturbed position is obtained by 
integration of the slope dq/ dx; this slope is related to v 2 by 

dx U 2 

Like u 2 , v 2 has been obtained as a Fourier integral involving G 2 (X) 
and Hg(X), and G 2 and H 2 are now known. Again, evaluation and 
simplification of the integrals is an involved process, and the final 
displacements are exhibited in series form and evaluated numerically. 


RESULTS AND DISCUSSION 

Present Modified Tsien-Finston Model 

Mach number par ame ter 9. - In the results herein the Mach numb ers 
M-j_ an d of the respective supersonic and subsonic streams appear in 
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a joint parameter 0 as well as individually. . For a given value of > 9, 
M 2 is a prescribed function of Mj (see equation (B32)). TJie curves of 
M 2 against Mj_ for several values of 9 are given in figure 2, which 
has been reproduced from figure 3 of reference 2. Tsien and Finston 
note that if the representative Mach numbers for the boundary-layer 
flow and the free stream are taken to be 0.8 and 2, respectively, then 
the corresponding value of 0 in the flow model is approximately 3it/4. 

Pressure distribution . - Curves of the pressure distribution along 
the upper interface calculated from equations (B62) and (B63) with the 
aid of functions tabulated in reference 2 are given in figures 3(a), 

3(b), and 3(c) for values it/4, jt/2, and 3rt/4, respectively, for the 
Mach number parameter 0. These curves exhibit the same infinite pres- 
sure peak at the point of shock incidence (£ = 0) found by Tsien a nd 
Finston with the original model. A similar peak was found by Lighthill 
(reference 5) in his model with continuous, rather than stepwise, veloc- 
ity profile; in that case the peak occurred inside the simulated boundary 
layer at the point of reflection of the shock from the sonic line. Such 
a peak will be modified in a real flow and probably will appear as a 
localized region of high pressure. 


Displacement of upper and lower interfaces . - The respective dis- 
placements of the upper and lower interfaces from their undisturbed 
positions have been designated and q L ; both are functions of the 

longitudinal distance x. Convenient corresponding nondimens ional 


parameters are 


% 


2bp 2 ' 2t,’ 2bp 2 2e 


and 


Jtx 

2bp 2 ' 


Curves of the displacement 


against longitudinal distance in terms of these nondimens ional parameters 
are plotted in figures 4(a), 4(b), and 4(c) for 0 equal to */4, it/2, 
and . 3it/4, respectively. The infinite slope at x = 0, associated with 
the pressure peak there, is so localized on the drawings as to be 
scarcely distinguishable. 


The scales are such that the vertical and horizontal displacements 
of a given interface are to the same geometric scale for the case 
e = 1/8 radian; for this value of a the separation of the interfaces 
is to correct scale when, in addition, p 2 = jt/20 for 0 = jt/4 (M-j_ = 2.5, 

Mg = 0.988); p 2 = jt/lO for 0 = */2 (M x = 2.7, Mg = 0.95); and p 2 = */6 

for 0 = 3jc/4 (M-j_ = 3.1, Mg = 0.85). Numerical values of the ordinates 
and abscissas are also tabulated in tables I and II. 


Comparison with supersonic flow above still air, no boundary layer . - 
The curves are now to be compared with the conjecture given in the intro- 
duction: namely, that in the simplified model of figure 1 the upper 

interface behaves, in the vicinity of the shock, substantially as if the 
entire boundary layer (represented by region 2) were replaced by a dead- 
air region (an extension of region 3); that is, that the upper interface 
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deflects through an angle 2e just as the dead-air interface does in 
sketch 3. The comparison shows that the curves in figure 4 approximate 
the conjectured behavior; both interfaces asymptotically approach the 
expansive turn through an angle 2 e. 

For each interface, the final asymptote does not intersect the 
initial asymptote (which is the 5-axis) at the point of shock incidence, 

5 = 0; instead the point of intersection lies somewhat downstream thereof. 
The point of intersection of the two asymptotes may be thought of as the 
effective "comer" for the turn through 2 6; its displacement amounts 
from one to several "reduced" boundary-layer thicknesses (reduced thick- 
ness = bfSg)* depending on the value of the Mach number parameter 0. 

The local deviation from the asymptotic behavior, that is, from the 
rectilinear comer turn of sketch 3, is l im ited essentially to several 
reduced boundary-layer thicknesses upstream and downstream of the effec- 
tive corner. 


Further Generalization of Tsien-Finston Model and Limitations 

The flow model of figure 1 is supposed to simulate only one aspect 
of shock boundary-layer interaction, namely, the local effect of the 
separated region in the vicinity of the point of shock incidence. An 
obvious further generalization, later found to be of questionable util- 
ity, is sketched in figure 5 . Here the dead-air region is limited to a 
finite length to approximate the "separation bubble," and on either side 
thereof the solid-wall condition used by Tsien and Fins ton is reapplied. 

When points S and R are both more than several reduced boundary- 
layer thicknesses from I, the point of shock incidence, their effects 
in the neighborhood of I would be expected to be small. This consider- 
ation serves to justify the use, for the neighborhood of I, of the model 
of figure 1 in which S and R are effectively removed to minus and plus 
infinity, respectively. 

There is another mentionable difference between the models of fig- 
ures 1 and 5: In figure 1 the perturbation pressure in the dead-air 

region is taken to be zero (for simplicity), whereas a positive pertur- 
bation pressure is intended in figure 5. The pressure rise in figure 5 
is effected by means of a weak shock or Mach wave incident just above S, 
resulting in the general upward slope of the boundary layer between S 
and I. This more realistic positive pressure perturbation in the dead- 
air region may be introduced into the model of figure 1 by imagining the 
weak shock associated with point S to occur at minus infinity to the 
left. Then the entire flow field will have an upward deflection (say e-^) 
corresponding to the flow deflection through this shock. As a conse- 
quence, the interface deflections of figure 4 will have superposed a 
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deflection t] = e^x 
as though each figure were rotated counterclockwise through the angle 
*1* 


.or 


3t 


2j(B' 


1 

2e 


JtX 

2rtBc 


2e h 


this is substantially 


As a generalization from a single aspect to a more complete picture 
of shock boundary -layer interaction, the model of figure 5 has serious 
shortcomings. The model by itself is insufficient to dete rmin e the posi- 
tions of separation and reattachment S and R; for this purpose the 
o> neglected viscosity will have to-be brought into the picture. Moreover, 

8 between I and R the actual, boundary-layer flow will probably become 

turbulent , and the neglected turbulent momentum exchange across both of 
the idealized interfaces will markedly influence the shape of the effec- 
tive interfaces and the location of R. The transition to turbulent 
flow will probably be decisive in determining R unless the incident 
shock is relatively weak, since otherwise the large abrupt pressure rise 
upon reattachment would be expected to reseparate a laminar layer but 
not a turbulent one. For these reasons it has not seemed worthwhile to 
make a serious effort to obtain an analytical solution for this model. 


Applicability of the Assumptions of Boundary-Layer Theory 

On the present model of shock boundary-layer interaction, the region 
in which the upper and lower interfaces show substantial curvatures and 
deviations from parallelism has been shown to be limited to one or two 
reduced boundary-layer thicknesses on either side of the point of shock 
incidence j a similar result is found for the original Tsien-Finston 
model. These results imply that it is only in this region that the com- 
ponent pressure gradient dp/dy is significant. Elsewhere the primary 
assumption of the Prandtl boundary-layer theory, namely, that dp/dy 
may be neglected, appears to be applicable. 

The situation seems to be as follows: Both models of shock boundary- 

layer interaction completely neglect viscosity, but include the inertial 
effects in both x- and y-directions, although in idealized fashion. The 
boundary-layer theory, on the other hand, includes viscosity and the 
inertia effects in the x-direction alone, the latter with idealization. 

The results of these models, despite their neglect of viscosity and 
other idealizations, suffice to delimit the small region directly under 
the incident shock in which the inertia effects in the y-direction must 
be included. It appears to be a proper inference that the Prandtl 
boundary-layer equations, laminar or turbulent as the case maybe, and 
with provision for interaction with the outer flow, may be applied 
upstream and downstream of this region, except possibly at the point 
of reattachment of the separated bubble . The separated bubble itself, 
except that portion directly under the shock, is included in the 
region of applicability. 
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Furthermore, the extension of the boundary -layer equations into the 
forbidden region, with the assumption of a sudden turn of the displace- 
ment surface through an angle 2e at the shock, appears to be justified 
in a practical analysis where accuracy in local flow details at the 
shock is secondary to the over-all picture; this conclusion is in sup- 
port of a phase of the treatments of references 6 and 7 . (if the 
incident shock is too weak to cause a separated bubble, the 2s turn 
at the shock must be replaced by another condition.) Present unknown 
elements in such an application appear to be the choice of the point of in ; 

transition to turbulent flow, if the boundary layer is not initially tur- S I 

bulent, and equations relating pressure distribution and displacement w 

thickness for a separated region of a turbulent or transitiQnal boundary 
layer. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 19, 1952. • 

[Note added in proof: In a newly published paper ("A Mixing Theory for 

the Interaction Between Dissipative Flows and Nearly Isentropic Streams," 
Jou l x. Aero. Sci., vol. 19, no. 10, Oct. 1952, pp. 649-676) Crocco and 
Lees have, in effect, supplied the missing relation between pressure dis- 
tribution and displacement thickness for a turbulent or transitional 
boundary’- layer by means of an assumption on the rate of mixing with the 
outer flow. In their work the boundary layer is replaced by a quasi-one- 
dimensional flow with the aid of quantities calculated by the Prandtl 
boundary- layer theory. There is no conflict between the affirmative con- 
clusions of the present paper on the applicability of the Prandtl theory 
and the position taken in the Crocco-Lees paper [] 
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APPENDIX .A 


SYMBOLS 

The following symbols are used in this report: 


Ajl(X), B 1 (x) 

e 2 (\), f 2 (x), 

®2\X ) } H 2 (X) 
3-1 > -^ 2 > * • • 


h 

f - 

6 

M 

P 

Ap 

R 

S = 
U 

u,v 
W = 

x,y 

z = 

Pi = 
p 2 = 
r 


Fourier coefficients appearing in the solutions for 
the supersonic stream 

Fourier coefficients appearing in the solutions for 
the subsonic stream 

Designation for nine integrals considered in the 
analysis 

thickness of subsonic region 

function of x and y representing incoming waves 

function of x and y representing outgoing waves 

Mach number 

static pressure 

change in static pressure 

reattachment point in figure 5 

BbXpgj also separation point in figure 5 

undisturbed velocity 

disturbance velocity components in x- and y-direct ions, 
respectively 

x + pj_y 

rectangular Cartesian coordinates 

x/r 




- 1 


-A - M£ 2 

parameter introduced in equation (B3) 
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e 


\ 

€ = 

1 ( 5 ) 


p 

* 


deflection angle of compression wave 
displacement of interface from undisturbed position 


Mach number parameter defined by cos 


6 

2 


variable of integration 

jcx/2bP2 

step function def in ed by l(£ ) 
density 

disturbance velocity potential 




K <o 
£= o 
K >o 


Subscripts : 

L value of given function evaluated at the lower interface, 

y = -b 

U value of given function evaluated at the upper inter- 

face, y = 0 

y = +0 l imi ting value of given function as y approaches 

zero through values greater than zero 

y = -0 limiting value of given function as y approaches 

zero through values less than zero 

1 denotes quantities in the supersonic region y > 0 

(excepting Ip) 

2 denotes quantities in the subsonic region -b < y < 0 

(excepting Ig) 

Superscripts : 


* 


indicates limit of given function as y -*■ Oj for 
example, f(y)* = lim f(v) 

r-*o 
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APPENDIX B 


ANALYSIS 

In the present report , the boundary layer has been replaced by a 
uniform subsonic stream of finite width b which is bounded on one side 
by fluid at rest and on the other side by a uniform supersonic stream 
of semi-infinite extent. For purposes of analysis, this configuration 
will be considered oriented with respect to a rectangular Cartesian 
coordinate system such that the positive x-axis is parallel to,, and in 
the direction of, the undisturbed flow and coincides with the undisturbed 
interfane between the subsonic and supersonic streams, and such that the 
origin denotes the point of incidence of the given shock wave with this 
interface (see fig. 1} . Thus, for -b <y < 0 the flow is subsonic, 
while for y >0 the flow is supersonic. Furthermore , a ll quantities 
in the supersonic region will be denoted by the subscript 1 and all , 
quantities in the subsonic region will be denoted by the subscript 2. 


Supersonic Region 


Governing equation and solution for Incoming waves . - In the super- 
sonic stream, the linearized differential equation for the disturbance 
velocity potential is given by 


Pi" 


d 2 ^ 


cbc^ 




= 0 


(Bl) 


and has as its solution 


= f (x + P x y) + g(x - p-jy) (B2) 

where f and g are arbitrary functions. The undisturbed flow has 
been taken in the positive x-direction so that waves represented by 
f (x + p x y) will be incoming waves while those represented by 
g(x - Ppy) will be outgoing waves emanating from the upper interface. 


For purposes of the analysis, the form of the function f will be 
specialized and taken to be 

fo , x + p^y < 0 


f(* + P x y) = < 


U x e 

•rPl 


-y(x + Pjy) 


4 


(B3) 


x + pjy > 0 




U-]_ is the undisturbed velocity in the supersonic stream anti 


where 

where y and e are parameters. If y tends to zero, the inc oming 
waves will degenerate into a single compression wave with deflection 
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angle e, and the point of incidence of this single compressive wave 
with, the upper interface will he the origin of the given rectangular 
Cartesian coordinate system. 


Now, by definition, the components u and v of the disturbance 
velocity are given in terms of the disturbance velocity potential 
by the relations 


u 


d® 

= 5 ? 



(B4) 


Therefore, the components of the disturbance velocity of the inc oming 
waves are given by 

f 0 , x + p-jy < 0 


u l, incoming < 


(B5) 


-Upe -r(x + Ppy), x + p^y > 0 

~ Q 6 

V. Pi 


and 


fo 


v l, incoming ' 


^ -U-,6 e 


-r(x + p,y) 


x + p-^ < 0 


x + p x y > 0 


(B6) 


Incoming and outgoing waves in Fourier-integral form . - It will 
prove useful to express the disturbance velocities of the incoming waves 
in the equivalent Fourier- integral form. To do this, it is noted that 
the components of the disturbance velocity of the incoming waves are 
proportional to the function 


h(w)= 



7 


y 


v < 0 


\T> 0 


(B7) 


where w = x + p^y, and that h(w) can be represented in an entirely 
equivalent form by 
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Furthermore, the potential due to the outgoing waves can also he expressed 
in Fourier -integral form, namely. 


g(x - Ppy) = Up 


sin X(x - p-j^y) + Bp(x) 



where Ap(X) and B 1 (x) are undetermined Fourier coefficients. Thus, 
after differentiation of the expression for . g(x - Ppy) to obtain the 

disturbance velocities due to the outgoing waves, the separate distur- 
bance velocities can be superposed and the Fourier- integral expressions 
for the components of the total disturbance velocities (of both incoming 
an ti outgoing waves) obtained. Evaluated at the upper interface, these 
disturbance velocity components will then be given by 


Un 



5 - COS Xx + -o a 

x 2 r + x 2 


sin Xx 1 dX + 


U n 


^ XA-^X) cos Xx - XBp(x) sin XxJ dX 


(BIO) 


and 


m 

r ( r - 

v i,y=+o * J [ r 2 +x 2 


T 

O cos Xx + 


X ^ \ 

5 - sin Xx ) 

+ x 2 J 


dX - 


Pl^L 



XAp(X) cos Xx - XB 1 (X) sinXxJ dX (Bll) 


Subsonic Region 


Governing equation and general solution in Fourier-integral form . - 
For the subsonic stream, -b < y< 0, the linearized differential equation 
for the disturbance velocity potential is given by 


2 

?2 




(B12) 
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and the solution. <p 2 can "be expressed, among many possible form s, as 




dX + 


dX 

(B13) 


■where E 2 (X), F 2 (X), & 2 (X), and H 2 (x) are also undetermined Fourier 

coefficients. Then, by retracing the steps just taken for the super- 
sonic case, there is obtained 



and 



|e 2 (A) cos Xx - F 2 (\) sin Xxjx dX + 

jorgOO COB Xx - H 2 (X) sin Xxjx dX 

(B14) 


jEg(X) sin Xx + F 2 (X) cos Xxjk dX + 

Jg 2 (X) sin Xx + H^X) cos Xxjk dX 

(B15) 


Boundary Conditions 


In these solutions for the subsonic and supersonic streams, six 
unknown Fourier coefficients Ap(x), B;j_(x), E 2 (x), F. 2 (x), Gg(x)., and 
H 2 (X) exist. However, these coefficients are determinable by boundary 
conditions at the upper and lower interfaces. The physical conditions 
which must be satisfied across these interfaces are: (l) equal static 

pressure, and (2) equal inclination of the flow. At the upper interface, 
condition (l) gives, in the small perturbation approximation. 


l 
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(i) -Pi u i u i,,y = 40 P2 U 2 u 2,y=-0 for 811 X 


(B16) 


At the lower interface, the perturbation pressure is taken to he zero; 
condition (l) therefore reduces to 


(II) U2,ya-1) = 0 ^ all x 


(B17) 


Condition (2) gives no information at the lower interface since the flow 
speed is zero below the interface; at the upper interface, this condition 
yields (small perturbation approximation) 


(III) v 2 ,y = _ 0 u x = U 2 for an 


(B18) 


Establishment of the four equations for the four Fourier coefficients . 
With use of equation (B14), boundary condition II implies 


U. 


P" 

, I |e 2 (X) cos ^- x ~ F 2 (X ) sin Xxjx dX = 0 

Jo J 


(B19) 


so that E 2 (X) and F 2 (X) vanish identically and q> 2 reduces to 
(p 2 = U 2 / sinh|xp 2 (y+b2j jjl 2 (X) sin Xx + H 2 (X) cos Xxj dX 

(B20) 


whence 


u 2 = u 2 


sinh ^ 2 (y + ' b )] [} G 200 cos Xx - XH 2 (x) sin Xxj dX 


(B21) 


and 



v 2 = u 2?2 coshjxp 2 (y+b)| Jxg^X) sin Xx + XH 2 (X) cos Xxj dX 


(B22) 


Substituting equations (BIO) and (B2l), the latter ev alua ted at 
y = - 0, into equation (B16), that is, boundary condition I, gives 
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2 I -e / l r . *■ 

Pl u J / “5 ? cos + "7 

x x |*Pi / Vr 2 + r + 


XA^CX) costa - XB^(X) sin Xx 


(XbPgjjxGgCx) cos Xx - XH 2 (X.) sin XxJ dX 


(B23) 


whence , upon equating corresponding Fourier coefficients and making use 
of the fact that 


p 2 U 2 2 M 2 2 

— “2 = T7"2 for Pi = p 2 
Pl u l M 1 


there is obtained 


(B24) 


Mp r v 1 

Ai(k) - G 2 (X) -2 sinh (Xbp 2 ) = — ^ 


(B25) 


%w - “2^5 ^2 (Xbs 2 ) = ^ ^jsr 

Sim.i3.arly, boundary condition III (equation (B18)) implies 


(B26) 


cosd (XbPg) sin Xx + XH 2 (X) cos Xxj dX 

° ^ 2 T cos Xx + sin Xxj dX - 003 aljl 


- r j cos Xx + 


so that 


Aid) + H 2 (X) cosh (Xbp 2 ) - jjjq- ^ ^ 2 
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and 


B X (X) - G 2 (X) 


$2 

Pi 


cosh (Xbp 2 ) = 


6 1 

*0i r 2 + x 2 


(B29) 


Hence , the problem is now reduced to one of evaluating the four 
undetermined Fourier coefficients ^ A]_(X), B^(X), G 2 (X), £ind H 2 (X) by 

means of the four equations (B25), (B26), (B28), and (B29) . 


Evaluation of Fourier Coefficients 


Solving equations (B25), (B26), (B28.), and (B29) simultaneously 
yields the following expressions for G 2 (X) and H 2 (/v): 


G 2(X) 


and 


/ 2 i 2 2 

sinh 2 (Xbp 2 ) + ^ cosh 2 (Xbp 2 ) 


Pi 

2 


5 


= ^1 * r 2 + X 2 *(m^) Sinl1 + 01 COSh 


(B30) 


h 2 (x) 


|"P2 

LPi 


2 2 , /M 2 2 \ 2 , 1 

cosh (Xbp 2 ) + — 7?J Bi rin (Xbp 2 ) 


-2e 


*01 ’ r 2 + X 2 


[P2 r 

[Pi x 


M 2 

cosh (Xbp 2 ) g sirih (Xb0 2 ) 


3 


(B31) 


In addition, if the Mach number parameter 9 is defined such that 8, 

M2 2 


9 

cos ^ = 


Mn 


'm 2 2 \ 2 p 2 2 


(B32) 


\M, 


Pi 


a Although simplicity is not necessarily thereby achieved, 9 is intro- 
duced so as to maintain a parallel presentation with that of reference 2. 
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then 


, 9 2 

0-) 


P 2 

sirih 2 (X.hp 2 ) + — g cosl12 = 

Pi 


/ 2K 2 

1 /m 2 2 \ 

2 Wy 


eosh( 2 XbP 2 ) - cos 0 


cos 


2 £ 


(B33) 


so that G 2 CM and 22 ( 1 .) are also given by 

,2 

I-. /Mo sf 

g 2 M 


and 


1 A'z s ) 2 

L 2 Wy 


cosh (2Xbp 2 ) - cos 0 


2 0 
cos 2 


) 


-2 6 


* p i * r 2 + I U \M X 2 


/V 


P2 


sinh (Xhp 2 ) + w— cosh (X.bp 2 ) 


(B34) 


1 fa 2 2 \ cosh (2XbP2) - cos 0 

H2(x) l?Wy ^ 


cos 1 


-2 6 1* [P2 

* p i ’ r 2 + x 2 |_Pi 


r , 


cosh (Xhp 2 ) 2 siDl1 (^^ 2 ) 

Mn c 


(B35) 


Similarly, hy means of the same equations, the following expressions for 
A-^(X) and B 1 (X) are also obtained: 


[ 2 / 

—g cosh 2 (Xbgg) + ^ — W sinh^ (XbPg )| 


e_ 

*P 


2 

- z~ "g y cosh 2 (Xbp 2 ) - y 

i r + x L?i \%7 x 


m 2 2 P 2 


sinh^ (Xbp 2 ) - 2 — g- — sinh (Xbp 2 ) cosh (Xbp 2 )j 


(B36) 
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and 


B x (\) jj^g) sinh 2 (Xhp 2 ) + ^2 cosh2 C^2)J 

p 2 

= 77^ H % II siDh (XbP2) cosh (XbP2) 



(^b02) + 


P2 2 

P? 


cosh 2 



(B37) 


Solution for Pressure Along Upper Interface 


With the Fourier coefficients Ap(X), „Bp(X), Gg(X)_, 

evaluated, the solution of the originally stated problem is now complete. 
However, there exist other quantities of interest such as the pressure 
distribution along the upper interface and also the displacement of each 
interface from its undisturbed position. These quantities will now be 
cons idered . 


If the change in pressure from the undisturbed value is denoted by 
Ap, then the linearized value is 


Ap = - pUu 


Hence, denoting quantities at the upper interface by the subscript U 
and those at the lower interface by the subscript L, equation (B2l) and 
boundary condition I give 


(A?)xj ■= 0 



slnh (XhPg) ^XGgCx) cosXx -Xpg(X) sin Xx| dX 


(B38) 


where GgCX) and HgOO are given by equations (B34) and (B35). Or, 
with use of equation (B24) and also since ' equation (B32) further implies 
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and 


$2 


sin 



(B40) 


equation (B38) can then he written as 

(Ap) T 


U 


I pA 2 



cos Xx + 


r 2 + x 2 


sin Xx) dX - 


y cos Xx + X sin Xx 

(y 2 + \ 2 /)(cosh 2Xbp 2 - cos 9 ) 


dX + 



X sinh (2Xbp 2 ) cos Xx - y sinh (2Xbp 2 ) sin Xx 
(r 2 + \ 2 )(cosh 2Xbp 2 - cos 9 ) 

(B41) 

Limit as y -» 0 (incoming waves -*■ step function) . - If T->0, the 
incoming waves will degenerate into a single compression wave with 
deflection angle t. Consequently, in order to obtain the value of the 
pressure change at the upper interface for a single compression wave of 

(Ap)y 

deflection angle e, let y-+0 in the expression for -n — . 

5 PA 2 


From equation (B8) a 


^ere, the subscripts under I designate the various integrals rather 
t han any position in a subsonic or supersonic stream. 
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x < 0 
x = 0 
x > 0 

(B42) 


x < 0 
x = 0 
x > 0 

(B43) 

For the other integrals of equation (B4l) , a theorem justifying the 
interchange of the processes of integration and of passing to a limit is 
necessary. One such appropriate theorem which will suffice for the 
analysis in this report will now be stated (reference 10, p. 474). 

DEFINITION . - Let R = (a® exp) denote the rectangle determined by 
a<x<«, a<t<p, p finite or infinite. The function f(x,t) is 
said to be REGULAR IN R when: (l) f (x,t) has no point of infinite 
discontinuity in Rj (2) f (x,t) is integrable in a < x < ® for each 
t in a < t < p. 

THEOREM I . - (l) Let f(x,t) be regular in R = (a® exp)/ p finite 
or infinite, except possibly on the lines x = a^, . . . , x = a^.. (2) Let 

f (x,t)dx 


be uniformly convergent in a < t <’p. (3) As t-*tQ, tQ finite or 

infinite, let f(x,t) approach a limiting function F(x) uniformly in 
any region a <x'<i), except possibly on the lines x = a-^, ... ., x = a^. 

Then, ~ . 

lim / f(x,t) dx exists. 

h-^O 




COB Xx + 


2 ..2 

r + x 


o 


sin Xx ) dX = < — 


ite 


-rx 


whence 



cos Xx + 





a 
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(4) Let F(x) "be integrable in any region a < x < b. 


Then, 



Thus, for 


*2 S 


X cos Xx 


I (y 2 + X 2 )(eosh 2bXpg - cos 0) 
first letting z = l/y so that 


dX (B44) 


*2 - 



cos (yxz) 


(z 2 + l)(cosh 2byzp 2 - cos 0) 


dz (B45) 


and then applying theorem I, there is obtained 


lim Io * -z- 

n * 1 - COS 0 

r-*o 


For 


'0 


dz 


jt 


Jt 


z 2 + l 


2(1 " cos 0) . . 2 £ 


, 9 v 

4 B±nr g 


(B46) 


*3 * 


X sin Xx 


0 Cr 2 + X 2 )(cosh 2bPgX - cos 0) 


dX 


(B47) 


and 


*4 S 


X s nnh (2bpoX) cos Xx 
q (r 2 + X 2 )(cosh 2bp 2 X - cos 0) 


dX 


(B48) 
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let 


S — 2bX^2 


and 


5 = 


JtX 


so that 


I 3 = 



S sin^ 

jt 


dS 


[r 2 + 


— ^ ^ ] (cosh S 

«>V/ 


- cos 0) 


S sl.nh S cos 


I 4 = 


it 


as 


'0 M + _B 


9 9 ) (cpsh. s - cos 0) 

4-b P2 / 


Then, 



sin 


SK 

jt 


S(cosh S - cos 0) 


dS 


and 



gg 

sinh S cos — — 

jt 

S(cosh S - cos 0) 


dS 


Lastly, 



X s i.nh ( 2 Xbft 2 ) 8 i 11 

(y^ + \ 2 )(cosh. 2 XbP 2 - cos 0) 


d\ = 


Solution for single incoming compressive wave . - Combini n g 
five limits yields 


(B49) 

(B50) 

(B51) 

(B52) 

(B53) 

(B54) 

0 

(B55) 

these 
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(Ap) 


u 


lim r 2 

r-° £ p x u i 


~ 4g 

Pi 


— 

cos 2 ^ . l(? ) 

- 1 

f 1 ) 

P 0 J 

_ 


\4 sin 2 g/ 


~ sin 2 0 

2rt 


sinM 

3t 


S(cosh S - cos 0 ) 



where l(£) is defined "by 


r 


o, 


i(5) = < i 


V 


1. 


£ < 0 
5 = o 
K >o 


(B57) 


Now, it can he shown hy contour integration that, for £ > 0, 

-(2n+i) g | e -(^)g e (^) ? 


i 


sin§£ 
Jl 


S(coBh S - cos 9) 


1 , 1 
- 0 2jt sin 0 
8 stn^ g 


2a + 2 - ~ 2a + 17 


e>o 


(B58) 


and 



sinh S cos 


K 


2it / S(eosh S - cos 0}^ 2at 


-(2n+l)5 


■f #)< 




2n + 2 - ^ 2n + ^ 


,5>o 

■^(B59) 


Ftirthermore, since the integr ands axe odd and even functions of £, 
respectively, it can imme diately he deduced that, for £ < 0, 


slaM . 

as = - _1 


2a / Swoosh S - cos 3 ) 
JO 


(2a+l)e 


8' sia‘ 


2 0 2s ala 0 




N 

e> ~ / _ e 

2n + 2 - - 2a + - 

JC 


, S< 0 


(B60) 
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and 


/ sinh S cos 

/ ^ « as 1 ) (2n+l)£ 


/ Swoosh S - cos 9) 00 St A_/ e 

2n + 2 - — 2n + — 

Jo n= '° 

jt 'rt 


£<0 (B6l) 


Therefore, with the substitution of these results into equation (B56), 
it follows that: 


lim 
r-*- 0 


(Ap)u 

1 TT 2 

2 p l U l 



K < o 


K > o 


(B62) 


(B63) 


Displacement of Upper Interface 


A quantity of primary interest is the shape of the distorted inter- 
face between the supersonic and the subsonic flows. The slope of this 
interface is related very simply to the already determined pressure at 
the interface, according to the following considerations. If the dis- 
placement from the undisturbed position is denoted by t^, this slope is 
given by 


dThj Vo 

dx“ s all 5 ,y=0 (B73) 

approximately. Now, applying the boundary condition of equal slopes bn 
both sides of the interface, equation (B18) gives 

v 2 V 1 

Un=u7 all 5, y=0 (B74) 


Furthermore, v^_ can be related to uj_. In the region £ < 0, only out- 
going waves exist, corresponding to the g-function in equation (B2); thus 
by differentiation of g. 
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v 1 = - 6^ £ < 0, y = 0 (B75a) 

In the region £ >0, there exists , in addition, the incoming weak shock 
providing a flow deflection t for x~* 0, corresponding to the f-function 
in equation (B2) ; this provides- an additional term such that 

lim v^_ = - 82. u 1 ~ 26U3_ S, > 0, y = 0 (B75b) 

r~+ o 


When these equations are combined, there results 

f U 1 
" U-l 


K < 0, y = 0 


y 

lim — = \ 
_ dx ' 

X-t- 0 


(B76) 


V 


P-L ir - 2e , ? > o, y = 0 


Then, by virtue of the relation Apg = - p-^U^u^, this may he written 

r 


„ d % y 

r-* o 



K < 0, y = 0 


(B77 ) 


- 2e, % > 0, y = O 


This is the desired relation between the local slope and the pressure of 
the interface. 

Upon substitution of the calculated pressure, equations (B62) and 
(B63), there results 


^ rk ( 2 
26 • a V 
— sm0 ^ 


2n+2 - -I? 


lim < 

• r-+o 


n=0 2n + - 


1 A 2n * D 5 


? < 0 


(B78) 



- 1 h 


K > o 


(B79) 
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Expressions for displacement of upper interface . - Denoting 
lim by - , equation (B50) gives 

r-*-o L " v 

dT iu* d %* d€ dll u* it 


dx d£ dx d£ 2 bp 2 
Thus, with the assumption that 

6 ) = 0 


(B80) 


(B81) 


the expressions for 


dT lU* 

”dx~ 


can be integrated and will yield 


it 


2bp 2 

while for £ > 0, since 

-0 

I 

0 ) = 

v 

there is obtained 


(Z; e) = % sin e £ - 




So ( 2n+2 . |) 2 ' 


5 < o 

(B82) 


<%*(?< 0; 0) 

<3? 


d£ + 



d Tu*te>0; 0) x 

gg d ? (BBS) 


V ( 5i 8) = |5 sin 8 ^ - ---TV 

n=0 fen+2 - 




- 2e 


s in 0 e 

xt=0 


-( 2n+ 1) 


fa * i) 


s 

2 


J+0 

(B84) 


Hence 


m 

2C Sl nflV 

1 

, - fa + f) 

1 e 

BlU C» / 

it l_y 


+ / s \ 2 " / 0\ 2 

n=0 

(2n+2 - -) 

(2n + — ) (i 2n + 


\ */ 

V v V *y 


- 2e£, ? 


> 0 < B85 > 
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Therefore , in summary. 



By inspection of equations (B82) and (B85), it is noted that 
%*(£? ®) is finite and continuous at £ = 0, while, on the other hand, 
equations (B78) and (B79) show that dr^/dx is infinite at £ = 0. 


Displacement of Lower Interface 


Finally, the expression for the displacement of the lower interface 
will he derived. If the displacement from the undisturbed position is 
denoted by q L , this displacement will have a slope given in linearized 
approximation by 


dT ^L _ v 2,y=-b 


dx U c 


(B86) 


or, with use of equation (B22), 


dTj 


_L 

dx 


/ XP2 ( g 2 (X) sin Xx + H 2 (X) cos XxjdX 


(B87 ) 


where- Gg(X) and H 2 (X) are given by equations (B34) and (B35). That is. 
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Limit as y-»0 (incoming wave step function) . - ConBideration will 
now "be given the various parts of the expression for dT^/dx as one again 

passes to the limit y = 0. For 


X 6 S 


T sinh (Xbp 2 ) sin Xx 


j 0 (r 2 + X 2 )(cosh 2bXp 2 - cos 0} 
first letting X = yz, there is obtained 


dX 


(B89) 


sinh (ybzpg) • sin (yxz) 


lim I g = lim . 

Y-+ 0 r~*0 / (l + z^)(cosh 2yzbp 2 - cos 0) 

0 


- dz = 0 (B90) 


Similarly, for 


Ir, = 


y cosh (Xbp 2 ) cos Xx 
1 ^ (y 2 + X^)( cos - tL 2XbPg - cos 0) 


dX 


(B91) 


lim Iy = lim 


y -> 0 y -*-0 



1 - cos 0 


cosh (yzbp 2 ) cos (yzx) 

(l + z^)(cosh 2yzbp 2 - cos 0) 


dz _ 1 it 

IT 1 - cos 0 * 2 


dz (B92) 



1 + z 
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so that 


For 


and 


lim I 7 - q 

y-> 0 4 sin^ 77 



X cosh (Xbf^) • sin 
(r 2 + X 2 )(cosh 2Xb6 2 - cos 0) 


dX 



X sihh (Xhp 2 ) . cos Xx 
(y 2 + X 2 )(cosh 2Xhp 2 - cos 0 ) 


dX 


introducing S and £ as before yields 


and 


whence 


and 




cosh - 77 sin 


it 


S(cosh S - cos 0) 


dS 



, ; s s? 

sinh 77 cos — r 

“ JC 

S(cosh S - cos 0) 


(B93) 


(B94) 


(B95) 


(B96) 


(B97) 


(B98) 


(B99) 


0 


2665 



2665 


NACA TN 2860 


31 


Furthermore, following reference 2, it can be shown that, for K > 0, 


r 


. a 2 e „ . „ i_, 

4 Bin ^ 2 sin g n.0 


1 ’ V ( ,'D (2n+l) 5 


-1) e x 




lim I 8 a ( 
r-*o 


* 1 V , ,*n -{2rn-l)E 

. : 2 e - ~~~s Lj e 

4 s1d £ ^ 2 sin ^ n=0 


2n + 2 - — 2n + — 

C.L- ^ ^ 




K< o 


2n + 2 - - 2n + - 
# jt 


(BXOO) 


£>o (bioi) 


and 


z'" 


1 E (-l) n e (2n+l)5 


2 cos ^ n=0 


.f?> -(¥> 


11m I. = < 
r-»o s 


2 cob 


(-1)“ e'^ 1 ^ 


n $ 

2n + 2 - — 2n + - 

_ jt it 


'-(^ e 


v. 


2d + 2 - £ 2n + — 

rt rt 


£ < o 


£ >0 


(B102) 


(B103) 


Simplification of expressions for slope . - Combining these limits 
and using equations (B 32) and. (B40) yield 


lim 
X~+ 0 


dT lL 

dx 


- 2e 



% >0 


(B104) 


That is. 


lim 

x~+o 


dx 


2c 



? > 0 

(B105) 
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Furthermore , ty s ymme try cons iderat ions of the various integrands , it 
follows that, for £ < 0 , 


d %, 

lim -r=- = - 2 s sin 


T -+-0 


dx 


i-ff E C-U n 


2 1 it sin ^ n =0 


.(# ;fi 3 ' 


P 0 

2 n +2 - - 2 n + — 
it it 


it sin g n =0 


E C-u n 


( 2 n+l)? 


M* M k 


2 n +2 - £ 2 n + £ 


it 


it 



That is. 


dr l 

lim = - 2 e 

r-*> ^ 


- Bta ! S ( ' 1)n 


( 2 


2th-2 


2 n + 2 - — 
it 


, £ < 0 (B107) 


Therefore, in summary, 

r 




( 2 

n e ' 


; 2 n +2 - - £ 


„ dT lL / 
lim - 3 — = \ 

r-o 


0 1 ’ 

2 n + 2 - - 
it , 


- 2 e (1 - - sin _ 

it o 


^~\ -kn + D? 

iE(-D n e . i 


n =0 




2 n + — 

it 


£ < 0 (B107 ) 


£ > 0 (B1Q5) 


g-inal expressions for displacement of lower interface . - Xf 
dm. dm* 

lim — *=■ is denoted "by — 5 — , it is noted that 
r -0 ^ 

d n L *. d "i, as d V jt 

dx d£ dx d£ 2 h 6 2 


(B108) 


so that integration of equations (B105) and (B107), along with the Boundary 
condition 


Tl L *(-~> 0) = O’ 


(B109) 



2665 


5C 


MCA TU 2860 


33 


yields, for £ < 0, 


557 lL*t« 


2bpg 


f k 


2n+2 - ® £ 


- 1) 


2~ > % <0 


(B110) 


For £ > 0, 


\~0 


ni*(5}0) = 


dr] T *fe<O;0) 
”di 


d? + 





dT] L *(?>O;0) 
dT 


<a? 


so that 


* v^; 0 ) - ^ sin 1 2 <-i) n 

n*0 


2bp 2 


-1 1 e 

^2- 


-( 2 n + |« 


(& + 2 - |J + |) (zn + g 


5 > o 


(Bill). 


Therefore, in summary 



(bho) 


5 > 0 (Bill) 


Both Ti L *(£j0) and diy^/dx are continuous at = 0. The con- 
tinuity of t] l *{ Q ) is readily seen -while the continuity of dT] L -*/dx: 

may he shown by letting x = 0/2rt in the formula for jt/sin juc given on 
page 208 of reference 10. 


o 
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TABLE I - PERTURBATION OF UPEER INTERFACE 


£ V . ^ _ itx 

2bp 2 26 against ^ “ 25p 2 


0 = rt/A> rt/2 , and 3it/4 


K 

ft=Jty^4 

&=rt/2 

6b=37t/4 

-5.0 

0.00001 

0.00008 

0.00028 

-4.0 

.00007 

.00035 

.00097 

-2.5 

.00093 

.00335 

.00632 

-1.75 

.00345 

.0103 

.01620 

-1.00 

.0132 

.0324 

.0416 

- .75 

.0208 

.0480 

.0584 

- .50 

.0335 

.0722 

.0820 

- .40 

.0410 

.0856 

.0944 

- .25 

.0562 

.1120 

.1180 

- .10 

.0804 

.1520 

.1517 

- .05 

.0925 

.1715 

.1675 

0 

.1108 

.1993 

.1891 

.01 

.1138 ' 

.2009 

.1856 

.05 

.1215 

.2018 

.1679 

.10 

.1253 

.1958 

.1408 

.25 

.1240 

.1564 

.0450 

.40 

.1169 

.1003 

-0.0620 

.50 

.1017 

.0566 

- .1377 

.75 

.0516 

-0.0686 

- .3363 

1.00 

-0.0137 

- .2117 

- .5450 

1.75 

- .2803 

- .7168 

-1.2120 

2.50 

- .6323 

-1.3002 

-1.9155 

4.0 

-1.5293 

-2.6075 

-3.3750 

5 .0 

-2.2368 

-3.540 

-4.3650 
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TABLE H - EERTUEBATIOU OF LOWER INTERFACE 


it 

2b6 2 2s 


against 


€ “ 2b£ 2 


for 


9 = it/4, it/2, and 3it/4 


£ 

Q=yjlj4= 

Q=lt/2 

a=3it/4 

-5.0 

-C 

.000013 

-0.00011 

-0.00073 

-4.0 

- 

.00007 

- .00050 

- .00254 

-2.5 

- 

.00100 

- .00472 

- .01653 

-1.75 

- 

.00368 

- .0144 

- .0420 

-1.00 

- 

.0134 

- .0436 

- .1057 

- .75 

- 

.0205 

- .0625 

- .1428 

- .50 

- 

.0309 

- .0890 

- .1918 

- .40 

- 

.0362 

- .1021 

- .2150 

- .25 

- 

.0459 

- .1250 

- .2550 

- .10 

- 

,0577 

- .1522 

- .3005 

- .05 

- 

.0622 

- .1623 

- .3173 

0 

- 

.0670 

- .1728 

- .3345 

.05 

- 

.0721 

- .1841 

- .3528 

.10 

- 

.0775 

- .1962 

- .2718 

.25 

- 

.0959 

- .2358 

- .4330 

.40 

- 

.1173 

- .2803 

- .5013 

.50 

- 

.1338 

- .3130 

- .5503 

.75 

- 

.1811 

- .4065 

- .6860 

1.00 

- 

.2389 

- .5168 

- .8385 

1.75 

- 

.4735 

- .9293 

-1.3803 

2.5 

- 

.7942 

-1.4453 

-2.0100 

4.0 

-1 

.6420 

-2.6725 

-3.4025 

5.0 

-2 

.3253 

-3.5775 

-4.4675 


2665 
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